A novel multi-variate immunological approach, reveals immune variation associated with environmental conditions, and co-infection in the koala (Phascolarctos cinereus).

External signs of disease are frequently used as indicators of disease susceptibility. However, immune profiling can be a more effective indicator to understand how host responses to infection may be shaped by host, pathogen and environmental factors. To better inform wildlife health assessment and research directions, we investigated the utility of a novel multivariate immunophenotyping approach examining innate and adaptive immune responses in differing climatic, pathogen co-infection and demographic contexts across two koala (Phascolarctos cinereus) populations in New South Wales: the Liverpool Plains (LP), and Southern Highlands to South-west Sydney (SHSWS). Relative to the comparatively healthy SHSWS, the LP had greater and more variable innate immune gene expression (IL-1ß, IL-6), and KoRV transcription. During extreme heat and drought, koalas from the LP displayed upregulation of a stress pathway gene and reduced adaptive immune genes expression, haematocrit and plasma protein, suggesting the possibility of environmental impacts through multiple pathways. In those koalas, KoRV transcription status, Chlamydia pecorum infection loads, and visible urogenital inflammation were not associated with immune variation, suggesting that immune markers were more sensitive indicators of real-time impacts than observed disease outcomes.

Expanding the known distribution of phascolartid gammaherpesvirus 1 in koalas to populations across Queensland and New South Wales.

Koala populations across the east coast of Australia are under threat of extinction with little known about the presence or distribution of a potential pathogen, phascolartid gammaherpesvirus 1 (PhaHV-1) across these threatened populations. Co-infections with PhaHV-1 and Chlamydia pecorum may be common and there is currently a limited understanding of the impact of these co-infections on koala health. To address these knowledge gaps, archived clinical and field-collected koala samples were examined by quantitative polymerase chain reaction to determine the distribution of PhaHV-1 in previously untested populations across New South Wales and Queensland. We detected PhaHV-1 in all regions surveyed with differences in detection rate between clinical samples from rescued koalas (26%) and field-collected samples from free-living koalas (8%). This may reflect increased viral shedding in koalas that have been admitted into care. We have corroborated previous work indicating greater detection of PhaHV-1 with increasing age in koalas and an association between PhaHV-1 and C. pecorum detection. Our work highlights the need for continued surveillance of PhaHV-1 in koala populations to inform management interventions, and targeted research to understand the pathogenesis of PhaHV-1 and determine the impact of infection and co-infection with C. pecorum.

Development of diagnostic and point of care assays for a gammaherpesvirus infecting koalas.

The recent listing of koala populations as endangered across much of their range has highlighted the need for better management interventions. Disease is a key threat to koala populations but currently there is no information across the threatened populations on the distribution or impact of a gammaherpesvirus, phascolarctid gammaherpesvirus 1 (PhaHV-1). PhaHV-1 is known to infect koalas in southern populations which are, at present, not threatened. Current testing for PhaHV-1 involves lengthy laboratory techniques that do not permit quantification of viral load. In order to better understand distribution, prevalence and impacts of PhaHV-1 infections across koala populations, diagnostic and rapid point of care tests are required. We have developed two novel assays, a qPCR assay and an isothermal assay, that will enable researchers, clinicians and wildlife managers to reliably and rapidly test for PhaHV-1 in koalas. The ability to rapidly diagnose and quantify viral load will aid quarantine practices, inform translocation management and guide research into the clinical significance and impacts of PhaHV-1 infection in koalas.

Future-proofing the koala: Synergising genomic and environmental data for effective species management.

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (~430 to 300 kya), a major climatic reorganisation that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of vulnerable populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g., guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

Differences in constitutive innate immunity between divergent Australian marsupials.

Understanding immunity in wildlife populations is important from both One Health and conservation perspectives. The constitutive innate immune system is the first line of defence against pathogens, and comparisons among taxa can test the impact of evolution and life history on immune function. We investigated serum bacterial killing ability (BKA) of five marsupial species that employ varying life history strategies, demonstrated to influence immunity in other vertebrates. The brushtail possum and eastern grey kangaroo had the greatest BKA, while ringtail possums and koalas had the least. These differences were independent of social structure, captivity status and phylogeny, but were associated with diet and body size. Sex and disease status had no effect on BKA in koalas, however potential for differences between wild and captive koalas warrants further investigation. The current study has provided a foundation for future investigations into how adaptive and innate immunity interact in marsupials from an eco-evolutionary perspective.

First evidence of deviation from Mendelian proportions in a conservation programme.

Classic Mendelian inheritance is the bedrock of population genetics and underpins pedigree-based management of animal populations. However, assumptions of Mendelian inheritance might not be upheld in conservation breeding programmes if early viability selection occurs, even when efforts are made to equalise genetic contributions of breeders. To test this possibility, we investigated deviations from Mendelian proportions in a captive metapopulation of the endangered Tasmanian devil. This marsupial population is ideal for addressing evolutionary questions in conservation due to its large size, range of enclosure types (varying in environmental conditions), good genomic resources (which aid interpretation), and the species' biology. Devil mothers give birth to more offspring than they can nurse in the pouch, providing the potential for intense viability selection amongst embryos. We used data from 140 known sire-dam-offspring triads to isolate within-family selection from population-level mechanisms (such as mate choice or inbreeding), and compared observed offspring genotypes at 123 targeted SNPs to neutral (i.e., Mendelian) expectations. We found lower offspring heterozygosity than expected, and subtle patterns that varied across a gradient of management intensity from zoo-like enclosures to semi-wild environments for some loci. Meiotic drive or maternal-foetal incompatibilities are consistent with our results, although we cannot statistically confirm these mechanisms. We found some evidence that maternal genotype affects annual litter size, suggesting that family-level patterns are driven by differential offspring mortality before birth or during early development. Our results show that deviations from Mendelian inheritance can occur in conservation programmes, despite best-practice management to prevent selection.

Characterization of reproductive gene diversity in the endangered Tasmanian devil.

Interindividual variation at genes known to play a role in reproduction may impact reproductive fitness. The Tasmanian devil is an endangered Australian marsupial with low genetic diversity. Recent work has shown concerning declines in productivity in both wild and captive populations over time. Understanding whether functional diversity exists at reproductive genes in the Tasmanian devil is a key first step in identifying genes that may influence productivity. We characterized single nucleotide polymorphisms (SNPs) at 214 genes involved in reproduction in 37 Tasmanian devils. Twenty genes contained nonsynonymous substitutions, with genes involved in embryogenesis, fertilization and hormonal regulation of reproduction displaying greater numbers of nonsynonymous SNPs than synonymous SNPs. Two genes, ADAMTS9 and NANOG, showed putative signatures of balancing selection indicating that natural selection is maintaining diversity at these genes despite the species exhibiting low overall levels of genetic diversity. We will use this information in future to examine the interplay between reproductive gene variation and reproductive fitness in Tasmanian devil populations.

A demonstration of conservation genomics for threatened species management.

As species extinction rates increase, genomics provides a powerful tool to support intensive management of threatened species. We use the Tasmanian devil (Sarcophilus harrisii) to demonstrate how conservation genomics can be implemented in threatened species management. We conducted whole genome sequencing (WGS) of 25 individuals from the captive breeding programme and reduced-representation sequencing (RRS) of 98 founders of the same programme. A subset of the WGS samples was also sequenced by RRS, allowing us to directly compare genome-wide heterozygosity with estimates from RRS data. We found good congruence in interindividual variation and gene-ontology classifications between the two data sets, indicating that our RRS data reflect the genome well. We also attempted genome-wide association studies with both data sets (regarding breeding success), but the genomic data suffered from small sample size, while the RRS data suffered from lack of precision, highlighting a key trade-off in the design of conservation genomic research. Nevertheless, we identified a number of candidate genes that may be associated with variation in breeding success. Individual heterozygosity, as measured by WGS or RRS, was not associated with breeding success in captivity but was negatively associated with litter sizes of breeding females in the RRS data set. Our findings enable conservation managers to have confidence in RRS data while understanding its limitations, and provide avenues for further investigation into which processes underlie variation in breeding success in captive Tasmanian devils. We caution, however, that deep functional insights using RRS may be impaired by a lack of precision, especially when marker density is low.

Too much of a good thing? Finding the most informative genetic data set to answer conservation questions.

Molecular markers are a useful tool allowing conservation and population managers to shed light on genetic processes affecting threatened populations. However, as technological advancements in molecular techniques continue to evolve, conservationists are frequently faced with new genetic markers, each with nuanced variation in their characteristics as well as advantages and disadvantages for informing various questions. We used a well-studied population of Tasmanian devils (Sarcophilus harrisii) from Maria Island, Tasmania, to illustrate the issues associated with combining multiple genetic data sets and to help answer a question posed by many population managers: which data set will provide the most precise and accurate estimates of the population processes we are trying to measure? We analysed individual heterozygosity (as internal relatedness, IR) of 96 individuals, calculated using four genetic marker types (putatively neutral microsatellites, major histocompatibility complex-linked microsatellites, reduced representation sequencing, and candidate region resequencing). We found no correlation in IR values across marker types, suggesting that various genetic markers reflect different aspects of genomic diversity. In addition, some marker types were more informative than others for conservation decision-making. Reduced representation sequencing provided the highest precision (lowest error) for estimating population-level genetic diversity, and most closely reflected genome-wide heterozygosity both theoretically and empirically. Within the conservation context, our results highlight important considerations when choosing a molecular technique for wildlife genetics.